Is the squeezing of relic gravitational waves produced by inflation detectable? 
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Grishchuk has shown that the stochastic background of gravitational waves produced by an inflationary phase 
in the early Universe has an unusual property: it is not a stationary Gaussian random process. Due to squeezing, 
the phases of the different waves are correlated in a deterministic way, arising from the process of parametric 
amplification that created them. The resulting random process is Gaussian but non-stationary. This provides a 
unique signature that could in principle distinguish a background created by inflation from stationary stochastic 
backgrounds created by other types of processes. We address the question: could this signature be observed with 
a gravitational wave detector? Sadly, the answer appears to be "no": an experiment which could distinguish the 
non-stationary behavior would have to last approximately the age of the Universe at the time of measurement. 
This rules out direct detection by ground and space based gravitational wave detectors, but not indirect detections 
via the electromagnetic Cosmic Microwave Background Radiation (CMBR). 
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I. INTRODUCTION 

The physical processes that took place early in the history 
of the Universe are well understood at characteristic times t > 
1 sec, where t denotes proper time after the big bang, in the 
rest frame of the cosmological fluid. However the processes 
that took place at times much earlier than this, in what is often 
called the "very early Universe" are not well understood [|l]]. 

A number of theoretical models of the very early Universe 
have been constructed which differ significantly from each 
other. In many of these models, phase transitions that take 
place as the Universe cools (or supercools) play an important 
role. While it can be difficult to calculate the observational 
properties of a given model, when this can be done their pa- 
rameters can often be chosen so that, at late times, the re- 
sulting cosmology is consistent with the known observational 
properties of the present-day Universe. A significant prob- 
lem is that there is often no way to distinguish between these 
different models based on strictly objective observational cri- 
teria. 

One of the most successful models of the very early Uni- 
verse is the so-called inflationary model. Among the many 
variations of this basic idea, the most generic and well-studied 
is the canonical model of "slow-roll inflation" [|l],||. In this 
model, the stress-energy tensor of the very early Universe is 
dominated by a vacuum energy term, which results in an ex- 
tended epoch of exponential (or near-exponential power-law) 
expansion. During this epoch, the visible part of the Uni- 
verse becomes extremely homogeneous and smooth, because 
the initial fluctuations in the energy density are redshifted into 
insignificance. 

Inflationary models have been extremely successful for sev- 



eral reasons. First, they solve two outstanding problems of 
modern cosmology, the horizon and flatness problems. Sec- 
ond, they can be described and understood with simple an- 
alytical models. Third, they make minimal assumptions. 
And finally, they make very definite observational predictions 
about the present-day properties of the Universe. In certain 
cases, for example in their predictions about the temperature 
anisotropics in the cosmic background radiation, these infla- 
tionary models are in excellent agreement with the observa- 
tional data [|| . However other models of the early Universe, 
such as the cosmic string model, are also in good agreement 
with this data [||. 

One of the ways in which different models of the very early 
Universe can (at least in principle) be distinguished is in the 
predictions that they make about the stochastic background of 
gravitational radiation [H,pp. This is a weak background of 
gravitational radiation, typically isotropic in its distribution, 
which is produced at very early times either by the large-scale 
motions of mass and energy in the cosmological fluid, or by 
the process of parametric amplification engendered by the ex- 
pansion of the Universe. Because gravitational forces couple 
so weakly, this radiation typically evolves free of disturbing 
influences, and at the present time, its properties provide a 
picture of the state of the Universe at very early times. To 
emphasize this point, consider that the present age of the Uni- 
verse is approximately 
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where the Hubble expansion rate today is 

Ho = lOOkmsec^^ Mpc"^ /iioo 



3.2 X lO-^^sec"^ h 



100- 
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Here ft-ioo ~ 0.65 is a dimensionless parameter. The cos- 
mic microwave background radiation (CMBR) provides us 
with a picture of the structure of the Universe at a time about 
t « 100, 000 years w 3 x 10^^ sec after the big bang. This 
should be compared to the present time Tq sa 3 x lO'^^ sec. 
On general grounds, the generation of ground-based gravita- 
tional wave detectors currently approaching completion might 
be able to detect a relic background of gravitational waves 
produced around t w 10^^^ sec after the big bang, providing 
us with a picture of the very early Universe [Bl). 

Of particular interest to us are the two LIGO detectors 
which will begin engineering shakedown within the next year 
[0], and the European VIRGO ||] and GEO-600 projects [|]. 
These ground-based detectors are sensitive in the frequency 
range from ^ 10^ Hz to ^ 10"^ Hz; long baseline detectors in 
space such as the proposed LISA [ |l(3| | project would extend 
the frequency range downwards by about four orders of mag- 
nitude to ^ lO^'^ Hz. The exciting prospect is to use these 
instruments to learn something about the physical processes 
that took place in the very early Universe. 

Models of the very early Universe typically predict a relic 
background of gravitational waves which is stochastic, in the 
sense that the gravitational strain at any point in the Universe 
is a non-deterministic function of time which can only be char- 
acterized in a probabilistic way. Often, the gravitational strain 
arises from the sum of a large number of independent pro- 
cesses and so the central limit theorem implies that the result- 
ing random process is Gaussian. In addition, in many models 
the stochastic background is stationary to a good approxima- 
tion. 

Grishchuk has shown that the stochastic gravitational wave 
background produced by inflation is Gaussian, but is not sta- 
tionary Jll|-13|]. This is because, in inflationary models, 
the gravitational waves are produced by a process of para- 
metric amplification of vacuum fluctuations, resulting in a 
squeezed quantum state today. While Grishchuk's derivation 
used the language and formalism of quantum optics, the non- 
stationarity has also been derived using the standard methods 
of curved space time quantum field theory 1 14,1^]. 

The non-stationarity of relic waves from inflation is signif- 
icant for two reasons. First, one calculation indicates that it 
might make the background easier to detect, because the in- 
tegrated signal-to-noise ratio can be made to rise faster than 
the square root of the observation time Jig], which is how 
it would rise for a stationary and Gaussian background [pT|]. 
Second, it provides in principle a definite test of inflation: if 
the squeezing is observed then it provides additional evidence 
in favor of inflation, and if the background is observed to be 
not squeezed, then it falsifies the theory. 

In this paper we present a detailed analysis of the detectabil- 
ity of the non-stationary statistical properties of relic gravita- 
tional waves. We show that, unfortunately, the effect can not 
be observed in an experiment whose duration is short com- 
pared to the age of the Universe at the time of the measure- 
ment, which rules out the possibility of measuring the non- 
stationarity with future ground-based and space-based gravi- 



tational wave detectors. The reason is a combination of two 
effects. First, the statistical properties of the background to- 
day vary over a frequency scale ~ T^^ 



10-1^ Hz, and 



any practical experiment (whose duration is of order a few 
years or less) will have a frequency resolution which is coarse 
compared to To ^. Second, there is an inherent and unavoid- 
able limitation in the accuracy of our measurement of statis- 
tical properties of a random process, due to the stochastic na- 
ture of the random process. This limitation is analogous to 
what is called "cosmic variance" in measurements of CMBR 
anisotropies. 

We also show that even in a gedanken experiment lasting 
long enough to distinguish between the stationary and non- 
stationary background, there does not appear to be any way to 
exploit the latter to make the integrated signal to noise grow 
faster than the square root of the integration time [ 16 1. 

We note in passing that it may be possible in the future 
to probe the squeezed nature of relic gravitons, albeit indi- 
rectly. Sakharov oscillations in the angular spectrum of tem- 
perature fluctuations in the CMBR are produced by a combi- 
nation of tensor perturbations (gravitational waves) and scalar 
perturbations. The existence of these Sakharov oscillations 
is directly related to the squeezed nature of the tensor and 
scalar perturbations [p^. If it is possible in the future to 
disentangle the scalar and tensor contributions to the CMBR 
anisotropies (perhaps using polarization information Jl^]), it 
may be possible to confirm that the gravitational waves are 
indeed squeezed. However it appears that will not be possi- 
ble with the upcoming MAP and PLANCK experiments Jl^]. 
Note that the possibility of such a measurement is consistent 
with our analysis, since CMBR photons effectively "measure" 
the gravitational perturbations over a time larger than or of or- 
der the age of the Universe at recombination. 

The paper is organized as follows. In Appendix ^ we re- 
view the predictions of general inflationary models for the sta- 
tistical properties of relic gravitational waves. In Sec. |lIA| we 
extract the cogent features of those predictions in the context 
of a simple model problem, a scalar field in 1 + 1 dimensions, 
and describe the differences between stationary and squeezed 
random processes in this context. In the remainder of Sec. || 
and in Appendices g]and y we analyze in detail measurements 
that take place at one point in space. We show that in order to 
distinguish between the stationary and squeezed processes in 
this context, a necessary but in general not sufficient condition 
is that individual modes need to be resolved with a precision 
that can only be achieved with an experiment lasting a sub- 
stantial fraction of the age of the Universe. The arguments 
do not rely on any specific inflationary model. Similar con- 
clusions were reached in work by Polarski and Starobinsky 
(page 389 of 0). 

In Sec. |l|we present a more detailed analysis that relaxes 
the assumptions of our simplified analysis of Sec. y. We con- 
sider an explicit model of slow -roll inflation, and calculate the 
two-point correlation function C(t, t') = {hi{t)h2{t')) of de- 
tector strain at two different sites. This is the quantity that 
will be measured in the upcoming generation of experiments 
in searches for a stochastic background [0,0]. We are able 
to obtain and confirm the main results of Ref. \uAn: due to 



the squeezing, non-stationary terms (not functions of t — t') 
arise in C{t,t'). Unfortunately we also show that the non- 
stationary terms are averaged away in any observation which 
is short compared to Tq. In addition, from the explicit ex- 
pression of the correlation function it can be seen that the in- 
tegrated signal to noise after optimum filtering should grow 
with the square root of the integration time [ p^ and not faster. 
Of course it is a dangerous game to claim that something is 
not possible. The history of physics is full of "no-go" theo- 
rems that have been circumvented by clever experiments. But 
our work does show that identifying the squeezed nature of 
the gravitational waves produced in an inflationary Universe 
is probably not possible with ground and space based gravita- 
tional wave detectors. 



II. SIMPLIFIED DETECTABILITY ANALYSIS 

A. Stationary and Squeezed random processes 

The statistical properties of relic gravitational waves, as 
predicted by inflationary models, are summarized in Ap- 
pendix ^ In this section we present an analysis of the de- 
tectability of the non-stationarity in the following simplified 
context. Consider a flat 1 + 1-dimensional spacetime with 
topology R X S^, with spatial circumference 2ttL, and metric 



ds^ 



L^{-dt^ + dx^ 



(2.1) 



where the periodic spatial coordinate x e [0, 2tt) and t is 
now a dimensionless time coordinate. We assume a compact 
spatial topology for technical convenience only; the discrete 
mode normalization is simpler than the corresponding contin- 
uum normalization. As a model of gravitational wave pertur- 
bations, consider solutions to the scalar massless wave equa- 
tion 
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(2.2) 



in the spacetime (2.1 ). The periodic nature of the spatial sec- 
tions means that solutions to the wave equation are periodic 
in timeQ and have discrete frequencies. The frequencies are 
separated by A/ = 1/{2ttL). 

The analog in this context of a stationary, Gaussian stochas- 
tic background (the naive prediction of inflationary models) is 
a scalar field of the form [ : 



* Modulo the uninteresting solutions which are linear functions of 
time and independent of x. We shall restrict attention to the periodic 
solutions, and consequently without loss of generality we can also 
restrict the range of the time coordinate to t G [0, 2tt) and imagine 
that the Universe has topology S^ x S^. 

^The n = terms have no gravitational wave analog and should be 
excluded from all formulae. 
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rn cos{\n\t + nx + (pn) ■ (2.3) 



This solution to the wave equation ( |2.2| ) describes a particular 
type of s tochastic random process. It can can be obtained from 
Eq. (A7) below by dropping the second term, by restricting 
the unit vector n to have two allowed values, and by using a 
discrete instead of a continuous mode normalization. Here the 
quantities r„ and(/)„ forn — ±1,±2, . . . are random variables 
whose statistical properties are given by 
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(2.4) 



where x„ and y„ are independent, zero-mean Gaussian ran- 
dom variables with (a;^ ) = {y^) = a'^/2. The variance cr„ 
can depend on n in an arbitrary way, but typically in infla- 
tionary models an has a power law dependence over a broad 
range of wavenumbers or frequencies; dn = a\n\^ for some 
constants a and /?. It follows from Eq. (2.4) that the r„'s 
and (pn's are all independent, that the r„'s are Rayleigh dis- 
tributed, and that the 0„'s are uniformly distributed over the 
interval [0, 27r). 

What Grishchuk has shown is that inflation leads to a 
slightly different random process, one which to very good ap- 
proximation can be written as 



'^squeezed (,''5 '^ ) 



E 



\/2r„ cos(ni) cos(ria; + 0„), (2.5) 



where the r„'s and 0„'s are distributed as before. The fac- 
tor of \/2 ensures that the processes ^-stationary and /igqucozcd 
have the same time-average d en ergy density. Note that the 
stationary random process (2.3) is invariant under changes 
t -^ i — Ai of the origin of the time coordinate, but that 
the squeezed random process (2^ ) is not. In Sec. |l| and Ap- 
pendix W we show that in inflation models, the simple form 
(2.5) is achieved when one chooses a particular origin for 
the confor mal ti me coordinate during the inflationary epoch. 
Again, Eq. (2.5) can be obtained from Eq. ( A7) below by re- 
stricting the unit vector n to have two allowed values, by us- 
ing a discrete instead of a continuous mode normalization, by 
taking x(fc) — 0, and by specializing the the limit of large 
squeezing |/3fe| > 1. 



Note that we can rewrite the squeezed random process ( 2.5 ) 



as 



/isqueezed(i, ^) = E ^ '"" cos(nt) cos(na; + (pn), (2.6) 



n=l 



where now the sum is only over positive values of n, and 
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It follows from Eqs. (2.4) and (2.7) that 
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(2.7) 



(2.8) 



where Xn and y„ are independent, zero-mean Gaussian ran- 
dom variables with (x^) — (j/^) = a^- Thus the squeezed 



random process can be described by half as many inde- 
pendent random variables as the stationary one. We also 
note that a solution to the two-dimensional wave equation 
is completely specified by h and h at some instant in time, 
where overdot denotes d/dt. Thus, given the two functions 

f{x) = /lstationary(0,a;) and g{x) = /lstationary(0,x) OnC Can 

uniquely determine r„ and </>„. The function /isqucozod has 
only half the number of degrees of freedom, and it is easy 
to show that it is a solution to the wave equation determined 
entirely ^ by f{x) = /isquoczcd(0, x). 

Below we shall be concerned with observations at a fixed 
point in space, which without loss of generality we can take 
to be the point x = 0. The two random processes evaluated at 
X — yield 

C30 

/istationary (i, 0) = ^ [Xn COs{nt) + {/„ Sm{nt)] 

71=1 

CO 

^squeezed (i, 0) = ^ \/2 S„ cos{nt), (2.9) 



where 



X X'fi -\- X—d 
Vn == -Vn + y-r 



(2.10) 



From Eqs. (2.7), ( |2.8| ) and (2.10) the random variables a;„, 
yn, in and y„ are all independent, zero mean Gaussians with 
variance cr,^ . 

Another convenient way of representing the processes 
/istationary(i, 0) and /isqucozod(i, 0) is to go from the time do- 
main representation of these fields to their frequency domain 
representation, through the Fourier transformation equations. 
The signals of Eqs. (E^ can be expressed as a superposition 
of a discrete infinite set of cosine functions: 

oc oo 

h{t, 0) = ^ i?„ cos(ni + 0n) = Re ^ i?„e'("*+*"\ 

n— 1 n— 1 

(2.11) 

defined, for every frequency n, by the Fourier amplitudes i?„ 
and phases 0„. These are what we shall refer to as "the Uni- 
verse modes". The relations between the coefficients appear- 



ing in Eq. (2.11 ) and those in Eq. (2.9) are, for the stationary 



case: 



Vn 



On = arg [xn + iyn] , (2.12) 

where the range of arg is [0, 2tt), and for the squeezed case: 

Rn ^ V2\Xn\ 



^It is determined up to an overall multiplicative constant, except on 
a set of measure of initial conditions. 



if S„ > 

. TT if Xn < 



(2.13) 



It is convenient to think of each mode as a vector in the 
complex plane having length i?„ and phase 6'„. This al- 
lows one to perform calculations using complex exponentials 
- which is easier - and also to have a useful pictorial repre- 
sentation of what the fi eld is. In the stationary case, it can 
be seen from Eq. (2.12) that the phases 0„ are uniformly dis- 



tributed on the interval [0, 27r). In the squeezed case, instead, 
the phases can only have two values, namely or tt. This 
means that whereas the mode vectors for the stationary case 
are randomly pointing in any direction in the complex plane, 
in the squeezed case they all lie on the real axis (see Fig. g). 
This is the origin of the name "squeezed". Also the distribu- 
tion of the amplitudes i?„ is quite different in the two cases, 
as can be seen from Fig. n^. 
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FIG. 1. Histogram of amplitudes R„ for the stationary (upper 
panel) and squeezed (lower panel) cases for a mode with ct^ = 1, 
from a total of 10"" points. 




stationary Squeezed 

FIG. 2. The distribution of the complex vectors of a mode with 
o"^ = 1 in the stationary and in the squeezed case (from a total of 
10^ points). 



B. Distinguishing between the two random processes in principle 



The question that must be addressed is this: what experi- 
ment can we perform that will in practice be able to distin- 
guish between these two different random processes, ( |2.3| ) and 
(2.5)? In other words, can we construct some observational 



statistic that will enable us to determine with high confidence 
which of these two random processes is present, given a single 
realization? 

To appreciate this issue a bit more clearly, consider that 
both /istationary and /isqucozod. ^Te solutions to the two- 



dimensional wave equation (p.2|). This fact can help us 
to compare the properties of these two random processes. 
Roughly speaking, the standard process /istationary is a sum 
of left and right moving waves with arbitrary amplitudes and 
phases: 
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(2.14) 



In contrast to this, the Grishchuk process /igqueezed is a sum of 
standing waves 



ft.squeezed(i, ^) = 51 ~7^ 1 COs[n(< + x) + (/)„] 
n=l V 2 I 

+ cos[7i(t — x) — 0„] > , (2.15) 



for which the left and right moving waves have equal ampli- 
tudes and correlated phases. 

Now suppose that we ask the following question: could we 
distinguish between these two possible forms of /i(i, x) based 
on unrestricted observations? The answer is obviously "yes". 
We would simply observe the evolution and see if there was a 
standing wave pattern or not: a standing-wave pattern would 
imply /isqucczod, anything else would falsify inflation. 

A slightly more difficult question is the following: sup- 
pose that instead we can only observe the function at a single 
point in space: /i(i, a; = 0). Again, could we distinguish be- 
tween these possibilities? Once again, the answer is clearly 
yes. By observing the process /i(t, 0) for < ^ < 27r and 
Fourier transforming, we can recover the complex coefficients 
i?„ exp[i6'„] in the expansion (2.11). If any of these coeffi- 



cients have non-zero imaginary parts, then the p rocess must 
be stationary rather than squeezed, by Eqs. ( 2.13 ). 

However, our assumption of spatial periodicity (which en- 
sures the periodicity in time) is just a computational device 
- essentially an infrared cutoff. We must at the end of our 
calculation let this cutoff go to infinity to obtain physical re- 
sults. In this limit, the measurement discussed above corre- 
sponds to a measurement of the random process over an infi- 
nite time, since the proper-time duration of the measurement 
is 27rL which diverges as L ^ cx). So, the conclusion so far 
is that it would be possible to distinguish between the station- 
ary and squeezed processes, provided that one had complete 



knowledge of the function h(t^x 
history of the Universe. 



0) over the entire time 



C. Distinguishing between the two random processes in practice 

Consider now realistic measurements of the stochastic 
background. Such measurements will be subject to three key 
constraints: 

• The function h(t^x = 0) can only be observed over 
a limited range of time Tg < t < Tg -V Tobs, where 
T ~ Ltis proper time, the starting time Ts is essentially 
the age of the Universe Tq, and the duration Tobs of 
the measurement is much smaller than Tq. In practice 
the longest possible observation times will be Tobs ~ 
1 year w 3 x 10^ sec. 

• The detectors available to us can not observe all of the 
modes, but just a (fairly narrow) range in frequency 
/ e [/o, ./o + A/] where typically /o ^ A/ ^ 100 Hz 
for ground based detectors, and /o ^ 10^* Hz, A/ ^ 
lO^'^ Hz for space based detectors. 

• The detectors that measure h{t, 0) are intrinsically 
noisy, limiting the accuracy of our measurements. 

We now show that the short observation time constraint 
(Tobs <C Tg ) alone makes it impossible to distinguish be- 
tween stationary and squeezed. In order to show this, we will 
derive the relation between the Universe modes and the ob- 
served modes that we actually measure. We will show that 
the latter are a weighted superposition of the former ones, and 
that the process of superposition makes it impossible to trace 
back the phase coherence of the squeezed case. 

Our analysis in this section will be limited to observations 
made at a single point in space. This assumption is of course 
not realistic, as realistic measurements will involve cross- 
correlating between spatially separated detectors. We present 
a more complete analysis which relaxes this restriction in Sec. 
|i| below. 

We start by making a convenient change in our representa- 
tion of the stochastic background. So far, we have treated the 
stochastic background as a continuous function of i - a ran- 
dom process h{t, 0) - whose frequency representation is dis- 
crete due to the assumed infrared cutoff (spatial periodicity). 
We now switch to a discrete representation of the stochastic 
background in the time domain, by modeling the background 
as a set of N numbers 



hj = h{jAt), 



(2.16) 



for < j < N — 1 where NAt = 2tt. Thus, we have effec- 
tively imposed an ultraviolet cutoff. At the end of our calcula- 
tion we can let At -^ and obtain a result that is independent 
of At [cf. Eq. (B20) below]. The discrete Fourier transform 



(DFT) of the sequence (2.16) is given by the relation 
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fe=0 



(2.17) 



whose inverse is 

hk 



N-l 
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(2.18) 



j=o 



We will also assume that the statistical properties of the quan- 
titie s hj a re the same as those of the quantities _R„ exp[i0„] of 
Eq. (2.11) above, namely 



{{Reh,y) = (l + e)<j 



{(Irah.Y) = [1 - e)aj 



(2.19) 



where e = in the stationary case and e = 1 in the squeezed 
case. Equation (2.19) will not be exactly true but will be true 



to an adequate approximation. 

Lets now take our measurement to consist of the quantities 



H, 



l-s+A 



(2.20) 



for < A < Af — 1, where s and M are fixed integers with 
< s < A^ and < M < N — s. Thus, the starting proper- 
time of the measurement is Tg = sLAt and the duration of 
the measurement is Tobs = MLAt. Henceforth capital Ro- 
man indices (A, B, C, . . .) will be measurement indices that 
run over 0,1, . . . ,M — 1, while lower case Roman indices 
(i, j, k, . . .) will run over 0, 1, . . . , A^ — 1. The DFT 
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(2.21) 



of the quantities (2.20) define the amplitudes of what we will 



call the measured modes. 

We can c ompu te the relation between the measur ed m ode 
amplitudes ( 2.21 ) and the Unive rse mode amplitudes ( 2.17 ) by 
inserting Eqs. (2.18) and (2.20) into Eq. (2.21), which gives 



Here the coefficients Wap are given by 



(2.22) 
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(2.23) 
(2.24) 



Equations (2.22) and (2.24) say that the complex vectors 
Ha (amplitudes of the observed modes) that define the DFT 
from the shorter time baseline ( ^.20| ) are a weighted sum of 
the vectors hp (amplitudes of the Univers e mo des) that define 



the DFT from the longer time baseline (2.16). The weights 
are complex numbers whose magnitudes depend only on the 
distance in frequency between A and p. In order to make this 
statement clearer, let us take a closer look at Wap and intro- 
duce a new index variable 



A'-A^ 



(2.25) 



This variable can take values from to M — 1, as does the 
index A, but is not an integer It measures the position of 
the points p on the A axis. In terms of this new variable the 
weights take the form 

TiA _ J_ -2mA's/M j7T(A-A')(i-i/Af) sin [7r(A — A )] 
^^'"^ s^ifliA-A')] 

(2.26) 

Note that when A' is an integer, we have Waa' = ^a'A^ 
which is necessary since the random sequences hj and Ha 
are real. Also we can write 



WAA'^^!M{\A-J^\)e'^--' 



with 
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siuTra; 
M sin jjx ' 



(2.27) 



(2.28) 



which shows that the absolute values of the weights depend 
only on the distance [A — A') between A and p. Figure 
shows the magnitudes JM^aA'! of the weights for the case 
M = 10, iV = 100. This figure makes it clear that the 
Universe modes tha t con tribute the most to a given observed 
mode, in the sum (2.22), are those which are nearest to it. 



It also shows that for any measured mode A, the amplitudes 
of the Universe modes to the right of A (i.e. modes with 
A' > A) are weighted in the same way as the modes to its left 
(A' < A), which follows from JM^a a+aa| — \Waa-aa\- 




FIG. 3. This plot shows the behavior of ] W^a' I as a function of 
A ~ A' for the case M = 10, A^ = 500. The maximum value is 



M 

AT 



10 
500' 



The crucial property of the weighting factors Waa' is that 
they contain oscillations on two different frequency scales. 
First, there is above-discussed oscillation in the magnitude 
IW^AA'I with characteristic frequency scale ~ T^].. Second, 
there is in addition a variation in the phase of Waa' , encoded 



in the factor exp[~2TnA' s / M] in Eq. (2.26), with character- 
istic frequency scale ~ T~^. 

We can now give the key argument. Consider a particular 
measured mode amplitude Ha, which corresponds to a phys- 
ical frequency /a — A/( MLA t). This mode amplitude re- 
sults from a vectorial sum ( 2.22 ) of Universe mode amplitudes 



hj. The Universe modes hj^pjA'/M which contribute signif- 
icantly to Ha all lie with a frequency interval around Ja of 
width r^ ^obs' because of the oscillations in JWatI'I- Now, 
if the weighting factors Waa' were all real in this frequency 
interval, then since in the squeezed case the Universe mode 
amplitudes hj are all real, then the resulting measured mode 



amplitude Ha obtained from the sum ( [2.22| ) would also be 
real. Hence, the observed mode amplitude would share with 
the original Universe mode amplitudes that property that its 
phase is constrained to be or n, and so the squeezing would 
be easily observable. 

Thus, we see that the crucial point is the extent to which the 
weighting factors can be treated as real in the frequency inter- 
val of width ~ T~^^ around the central frequency /a ■ Con- 
sider first the short-measurement regime Tobs <C Ts. In this 
regime, the phase of WaA' winds around between and 27r 
roughly Tg/Tohs ^ 1 times in the fr eque ncy interval, due to 
the factor exp[~2TTiA' s / AI] in Eq. ( 2.26 ), and so cannot be 
treated as real. In this case the observed mode Ha has a phase 
that is very nearly uniformly distributed over [0, 27r); the vec- 
torial sum (2.22) almost completely erases the peculiar phase 
behavior of the non-stationary process. This is in part due to 
the fact that the vectorial sum entangles phase and amplitude 
information and that the amplitudes of the Universe modes 
(although not the phases) are stochastic. In Appendix be- 
low we demonstrate this by deriving the distributions of the 
phase and amplitude of the observed amplitude Ha- 

So far, we have argued that the effect of the squeezing on 
the statistical properties of each observed mode Ha is very 
small in the regime Tobs ^ Tg. This gives the essential reason 
why the squeezing is not observable. However, we still need 
to close two loopholes in the argument. First, we have not 
yet shown that the effect of squeezing is unobservable, as we 
have not addressed the question of how accurately the statis- 
tical properties of each mode can be measured. In the case of 
the CMBR, for example, the anisotropics are smaller than the 
homogeneous background CMBR by a factor of 10^, yet those 
anisotropics are still observable and carry a great deal of infor- 
mation. We need to show that the intrinsic limitations in the 
possible accuracy of measurement of the statistical properties 
of each mode ("cosmic variance") are larger than the effect of 
squeezing on each mode. Second, it is insufficient to perform 
an analysis that focuses on each individual measured mode 
Ha one at a time, as the individual measured modes are not 
statistically independent, and we need to show that squeezing 
cannot be observed in any measurement that combines infor- 
mation from all the measured modes [ 20 1 . These two concerns 
are addressed and resolved in Appendix ^, where we demon- 
strate that the effect of squeezing is indeed not observable. 

For completeness, we now consider the long measure- 
ment regime Tobs ^ Tg, although as explained in the In- 



troduction this regime cannot be realized in practice. In this 
regime, the oscillation in the ph ase of Waa' due to the factor 
exp[—2TTiA's/M] in Eq. (2.26) is negligible. There is a phase 
oscill ation due to the factor exp[Tn{A — A') (1 — 1/^/)] in Eq. 
(2.26), which gives a total change of phase of order unity over 
the frequency interval. The resulting observed mode ampli- 
tude Ha is thus not constrained to be real. It turns out that the 
probability distribution for its phase ArgHA is exactly uni- 
form in the case when the spectrum is white, i.e., when <tj — 
const. When the spectrum is colored, the detectability of the 
squeezing depends on how large Tg is compared to a correla- 
tion time determined by the spectrum ctj (see Appendix g|). 



III. TWO-DETECTOR CORRELATION EXPERIMENTS 
A. Introduction 

In this section we extend the analysis of Sec. pi of the de- 
tectability of the non-stationarity to incorporate a number of 
complicating effects, including the effect of having measure- 
ments at more than one point in space. We restrict atten- 
tion in this section to the realistic case of short observations, 
Tobs ^ Tg. 

The standard method which will be used to detect a stochas- 
tic gravitational-wave background is based on correlating two 
widely-separated detectors. In the standard treatments of this 
technique, one considers a stochastic background which is sta- 
tionary and Gaussian. In this case it is easy to show that the 
integrated two-detector correlation arising from the stochastic 
background is proportional to the integration time, whereas 
the terms arising from the (assumed uncorrected) detector 
noise are proportional to the square root [n6p of the integra- 
tion time. 

As stressed by Grishchuk, the stochastic background pro- 
duced by an inflationary epoch in the early Universe is in a 
squeezed quantum state produced by the period of rapid ex- 
pansion. The statistical properties of this state are Gaussian, 
but not stationary. Unfortunately the non-stationary behavior 
does not make the stochastic background produced by infla- 
tion easier to observe than a similar backgrou nd wh ich is sta- 
tionary and Gaussianpj We show this in Sec. 



IIIB 



below by 

deriving an expression for the correlation function for the de- 
tector strains at two separated sites. We show that, in contrast 
with the correlation function for a stationary Gaussian back- 
ground (which depends only upon t — t'), the squeezed na- 
ture of the inflationary background gives rise to an extra term 



^Recent work by Turner [gjj has shown that within the context 
of slow-roll inflationary models, current observational limits on the 
CMBR imply that the gravitational wave stochastic background of 
inflationary models would be too weak to detect with either the cur- 
rent LIGO-I or planned upgrade LIGO-II detector, and is just below 
the limits of sensitivity of the proposed LISA [|lO| experiment. 



in the correlation function. This extra term explicitly mani- 
fests the non-stationary nature of the squeezed state: it is not 
a function of t — t' alone. Thus Grishchuk is entirely correct 
that the squeezed quantum state has different properties than a 
stationary and Gaussian background. However we show that 
observing this extra non-stationary term requires a frequency 
resolution of order 10~^^ Hz. Such incredible frequency res- 
olution is only obtainable in a gedanken experiment lasting as 
long as the present-day age of the Universe. 

The calculation in this Section makes frequent use of re- 
sults from two references. The first reference examines the 
gravitational waves produced by an inflationary epoch in the 
early Universe [p2|], and calculates the correlation function of 
the temperature anisotropics that they induce in the CBR. The 
frequency range of the waves influencing the observable part 
of the CBR anisotropy is between 10"^* Hz and 10"^'' Hz 
today. But the methods of Ref. ^2j\ can also be applied (as 
we do here) to calculate the properties of the waves at much 
higher frequencies: the lO^'^ Hz relevant to space-based or 
the 10^ Hz relevant to earth-based gravitational-wave detec- 
tors. The second reference that we make use of is a review 
article on how correlation techniques may be used to detect 
a stochastic background ^]. A reader wishing to follow our 
calculations in detail is advised to have copies of these two 
papers in hand. 

The calculation done in this Section is for an explicit model 
of slow-roll inflation, in the limit that the resulting spectrum 
is "flat" and not tilted. After completing the explicit calcula- 
tion of the correlation function, it is straightforward to show 
that in more general inflationary models, which produce a 
"tilted" spectrum, the same general conclusion holds: the non- 
stationary behavior can not be observed in any practical exper- 
iment. 



B. Calculation of the correlation function C{t, t') 

The most direct way to understand how the squeezing pro- 
duced in an inflationary model affects the correlation between 
separated detectors is to examine the two-point correlation 
function for the gravitational wave strain h at two detector 
sites. This is defined by: 



the strain at site i is given by: 



1 



C{t,t')^{hi{t)h2{t')) 



(3.1) 



where the subscripts 1 and 2 refer to the detector sites, and t 
and t' denote the time at these sites. The angle brackets can 
be given two possible (equivalent) meanings. They can de- 
note the average over a statistical ensemble of many different 
inflationary Universes, each starting from somewhat different 
initial conditions, but made statistically similar by a long pe- 
riod of exponential inflation. Equivalently, they can denote the 
expectation value in a quantum state that is a good approxima- 
tion to the present-day state of the Universe. 

The detector strains are given in terms of the metric pertur- 
bations of a Friedman-Robertson- Walker (FRW) cosmologi- 
cal model. Denoting the space-time metric by 



h{t) = -hab{t.^^){XtXl 



YtY^) 



(3.3) 



where x^ is the spatial location of the i-th site and Xf and Y°- 
are unit-length spatial vectors (with respect to the metric 8ah) 
along the directions of the orthogonal detector arms at site i. 
Our model of the Universe is specified via the cosmological 
scale factor a(j]) given as a function of conformal time r\ by 
Eq. (3.7). Note that we frequently use cosmological time t 



rather than conformal time 77 as a coordinate; they are related 
by dt = a{rj)dri. 

To compute the correlation function, we need an expression 
for the metric perturbation hab- This quantity is not deter- 
ministic. In classical calculations, it may be treated as a real 
stochastic random variable with certain statistical properties. 
For inflationary cosmological models, a slightly different ap- 
proach is useful, in which hab is replaced by the field operator 
hab of the linearized gravitational field. To see why this is 
both simple and useful, it is helpful to briefly review the fun- 
damental mechanism through which inflationary models erase 
the effects of the conditions present before the inflation begins 
and give rise to a spatially flat and homogeneous Universe to- 
day. 

In inflationary models, the Universe undergoes a period of 
(nearly) exponential expansion of the scale factor a. During 
this period, the energy-density of any matter or radiation is 
red-shifted away exponentially quickly. For massless parti- 
cles, the energy of a particle is proportional to a^^ and hence 
the energy-density in massless particles redshifts as a^*. For 
massive particles, any kinetic energy is quickly red-shifted 
away, leaving only the rest-mass behind. This energy is di- 
luted by the expansion of the spatial three-volume, so that 
the energy-density of massive particles redshifts away as a^^. 
What is left behind after the period of exponential inflation 
is just the vacuum energy associated with the cosmological 
constant; it is the dominant form of energy at the end of an 
inflationary epoch. 

The perturbations away from absolute uniformity in infla- 
tionary models are dominated by the effects of the zero-point 
vacuum fluctuations of the fields at the start of inflation [|23|]. 
Twenty five years ago Grishchuk showed that during the sub- 
sequent expansion of the Universe, the quantum zero-point 
fluctuations of the gravitational field modes are parametrically 
amplified [|4[|. This prediction has been subsequently con- 
firmed in many different ways, and the process of parametric 
amplification, which applies to both scalar and tensor pertur- 
bations, is one explanation of how the early Universe can gen- 
erate a spectrum of perturbations. 



ds^ = a^i-q) {-duf + [Sab + hab)dx°'dJ') 



(3.2) 



The above discussion, combined with Eqs. ( [3.1| ) and ([3.3|), 
implies that the correlation function today is well approxi- 
mated by the vacuum expectation value 

C(i,t') = -(0|/la6(t,Xi)/lcd(i',X2) + /lcd(t',X2)/lafc(t,Xi)|0) 

{X1X\ - Y{^Y,'){XIX^ - %%'), (3.4) 

where |0) is the (unique Hadamard de Sitter-invariant) quan- 
tum vacuum state of the inflationary Universe [ESI], and hab 



is the Hermitian field operator of the linearized gravitational 
field. 



The correlation function ( 3.4 ) can be calculated exactly, but 
since space and ground based detectors will only be sensitive 
to gravitational waves in the frequency range / G [lO^'^, 10"^] 
Hz, a number of simplifying approximations can be made. 
We do this by approximating the field operator within this 
frequency range. In general, the field operator is given by 
Ref. [^ as 



(3.5) 



habiv,^) = / d^k (e'^''[eab{k)(f>R{T],k)aR{k) 



+eli,{k)^Liv,k)dLik)] +h.c. 



where h.c. means Hermitian conjugate. This is a 
sum over wavevectors (labeled by k) of left- and right- 
circular polarization tensors e^h and 6*^^ and wavefunctions 
4'l.r{V7 k) exp(ik • x). To approximate the field operator we 
need to relate present-day frequency / to fc = |k| = \/kak°-. 
Since k always appears in the form Gxp(ik • x), from the form 
of the metric (3.2) one can see that the frequency / today is 
related to k by 



/ 



k 



2na{rio 



(3.6) 



We will now express this relationship between / and k in 
terms of the present-day Hubble expansion rate Hq, which 
will yield Eq. (3.11). This leads to relations (3.13) needed to 
approximate the field operator ( 3^ ) and to evaluate the two- 
point correlation function (3.4) in the frequency range of in- 
terest. 

The inflationary cosmological model is defined by the scale 
factor airj) in Eq. (4.1) of Ref. [|2|]: 



(2";^)"V^i) 



-oo < rj < rji de Sitter, 



a{r]) = ^ ij^O'iVi) m <'n <V2 radiation, 

. i(l + ;^)^m°(''i) ^2 < ?7 matter. 

(3.7) 

This scale factor describes three epochs: a de Sitter infla- 
tionary phase, followed by a radiation-dominated and then a 
matter-dominated phase. Letting t be the cosmological time 
defined by dt — a(rj)drj, the Hubble constant today is given 
by: 



Hn 



nda\ 



1 da\ 



'^V today U^c^Ty^^^^y 
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m + m 



{■>1Q 



(3.8) 



Here ?7o is the present-day value of the conformal time and 
r}2 is the value of the conformal time at the beginning of the 
matter-dominated phase of expansion. The redshift Z^q at 
which the matter and radiation energy densities are equal is 



a{m) 4 



10^ 



(3.9) 



Hence 770 is about two orde rs of magnitude greater than 772 and 
the Hubble constant today (3.8) is well approximated by 



i^o 



'70 



-0-^(770) « 10- 



-18 



Hz. 



(3.10) 



Having specified the cosmological model, one may now ap- 
proximate the field operator hab- 

We will be making approximations valid for the range of 
frequencies that might be obse rved by earth- or space-based 
detectors. It follows fromEqs. (3.6) and (3.10) that k is related 
to the present-day frequency / by 



47r/ 
'70^0' 



(3.11) 



In the frequency range of interest, krjQ w Airf/Ho ^ 10^^ 
is much larger than unity: krjo >> 1. From examination of 



Eq. ( p.9| ) this implies that one has kri2 ^ 10^^ so kri2 >> 1. 
Finally, from Eqs. (4.2) and (4.3) of Ref. @, if the period 
of inflation creates sufficient cosmological expansion to solve 
the horizon and flatness problems, then 



m 
m 



i + z, 



1 



^<io 
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(3.12) 



Here Zend ^ 10^^ is the redshift at the end of the inflation- 
ary epoch. Hence krji < 10^^ is much smaller than unity: 
h]i << 1. To summarize, we will approximate the field op- 
erator in the frequency regime where 



krjo » 1, 
h]i << 1, and 
k7]2 » 1. 



(3.13) 



A physical interpretation of these constraints will be given 
shortly. 

Any present-day detector correlation experiment will ob- 
serve the correlation function C{t, t') at the present time: t 
and t' are in the matter-dominated (prese nt-d ay) epoch. From 
the definition of the correlation function ( |3.4| ) and of the field 
operator ( |3.5| ) this means that we need an expression for the 
mode function (p{ri, k) in the present-day (matter-dominated) 
epoch, which is given by the final line of Eq. (4.17) of p^]: 



0(?7, k) = 0L,fl(?7, fc) = a ^l^aU'?^ k) 



P4Llt{ri,k). 



The mode functions in the present-day matter-dominated 
phase are given by Eq. (4.7) of p^] 



'/•mal(»7,^) =0,nat (??: ^) 

V Svr pp k(ri + r]2) 



(3.14) 



in terms of a spherical Hankel function of the second kind 

(2) ^ 

h\ . The Planck density is denoted by pp = ^^. 



The Bogoliubov coefficients a and /3 are given in terms of 
the corresponding coefficients for the transition between the 
de Sitter- and radiation-dominated phases and the transition 
between the radiation- and matter-dominated phases as 



a P 
13* a* 



a 



/3 
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(3.15) 



Since kr]\ << 1 and kri2 >> 1, one may make use of 
Eqs. (4. 12) and (4. 14) of Ref. ^ to obtain: 



-/3* 






(3.16) 
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Fina lly, m aking use of the definition of 0^^^^ given m 
Eq. (3.14), one obtains: 



• The gravitons are created in particle-antiparticle pairs. 
Since the antiparticle of a graviton is just a graviton of 
opposite momentum and helicity, gravitons are always 
created in oppositely-moving pairs 

• Since the amplitudes of these oppositely moving pairs 
of gravitons are exactly equal and their momenta are 
opposite [p7,28 1, they give rise to a pattern of standing 

waves. 

• This pattern of standing waves is apparent in the form 
of the mode function (j), whose complex phase is not a 
function of time (since (p is pure imaginary). 

These are precisely the conclusions reached by Grishchuk fol- 
lowing Eq. (7) of reference [pl|. 

The argument fc(ry — 771) of the cosine in Eq. (3.17) has a 



(j){r],k) Ki 4ik ^^^\ —miVo+m) ^ cos fc(?7 - ryi) 

V on pp 

4:ik^^''^\l -miln'^ cos ki] 

Sir Pp 



(3.17) 



where pds is the constant energy density during the inflation- 
ary de Sitter phase of expansion and pp is the Planck energy 
density. In deriving this expression we have made two approx- 
imations. During any period of observation lasting only a few 
years {v + V2) ^ (m + V2)- Also 



h^i\^) 



e t e 

X X : 



smce X 



Hv + ^2) ~ Hvo + V2) » 1 



Note that the field operator (3.5) is invariant under the com- 
bined transformations 



afl(k) -. e-'^'^'apik), 



,-iip(k)-^ 



4>{k,ri) 



jTp{k) 



(t){k,ri), 



(3.18) 



where ip{k) is an arbitrary function of k. The de Sitter vac- 
uum state is invariant under this transformation, since it is de- 
fined by a7<'(k)|0) = ai(k)|0) = 0. Hence without loss of 
generality we may multiply the mode function (3.17) by a k- 
dependent phase; the final physical results will not be affected 
by such a transformation. Changing this phase is analogous to 
multiplying the wavefunction of a quantum mechanical har- 
monic oscillator by a pure phase: it has no observable effects. 
In particular, changing the phase ijj is not equiva lent t o chang- 
ing in the argument of the cos appearing in Eq. ( 3.17 )^ 

The expressions f or the Bogoliubov coefficients ( 13.16| ) and 
the mode functions (3.17) have a number of interesting prop- 
erties: 

• The quantity |/3|2 is the (very large) number of quanta 
created by the "external" large-scale expansion of the 
Universe (or equivalently, by the parametric amplifica- 
tion of zero-point fluctuations) pq]. This is how in- 
flation gives rise to a potentially-observable stochastic 
background of gravitational waves today [ESh. 



simple physical interpretation. The number of cycles dN of 
a wave in the time interval dt at time t is dN — f{t)dt — 
f {rf)a{rf)drj = kdrj/2TT. This means that k{ri — rii)/2TT is 
the number of cycles of oscillation (in time) that the wave 
with wavenumber k has undergone since the end of the de 
Sitter phase at time rji . For frequencies observable by earth- 
and space-based detectors, the term krji << 1 and can be 
neglected. See also Appendix Q below. Later in this Section, 
we will encounter terms of the form k{r] zLrj'). 

Before completing the calculation of the two-point correla- 
tion function, it is useful to make a short digression. We will 
calculate the energy-density in gravitational waves using this 
formalism . Th e result illustrates precisely the effect predicted 
in Section II C : it is not possible to distinguish the the station- 
ary and squeezed states in local short-time observations. The 
energy-density in gravitational waves is given by 



Pgw ^ df 



dpgu 



df 

t^t' SzttG dt dt' 

ri^-q 6Z-nLi a[r]ja{r]'j drj drj 

The last term of the expression above, which is the two-point 
function of the field operator, may be deriv ed from the plane- 
wave expansion of the field operator (3.5 1 and the canonical 
commutation relations given by Eq. (2.17) of [G2]. Doing so 
yields 

{Q\hab{ri,^)hcd{ri' ,^m = j d^k e'^'^^-^'Uiv, kWiv', k) 
X [e,,(k)e:,(k) + e:,(k)e,,(k)] , (3.19) 

since the vacuum state |0) is annihilated by the operators dp 
and dL- Note that this quantity also appears in the definition 
of C(t, t') and will be referred to later in this section. Making 
use of the relationship between cosmological and conformal 
time dt — a{r])dr] and of the definitions (2.22-26) of [p2|], one 
obtains an expression for the energy-density in the stochastic 
gravitational wave background valid for the range of frequen- 
cies observable by earth- and space-based detectors: 
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k ^ cos^ fcr; dk. 



It is conventional to express this energy-density as a dimen- 
sionless spectral function ^gv/{f) which is the ratio of the 
energy-density in gravitational waves in a logarithmic fre- 
quency interval divided by the critical energy-density pc = 
■§;^Hq required to close the Universe^: 
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J O'Pgw "" i^as /-. 
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Note that ^gw{f) depends upon the time: if the statisti- 
cal properties of the state were stationary, it would be time- 
independent. For a single mode, the energy-density is an os- 
cillating function of time, exactly what we would expect for a 
standing wave, but without any spatial dependence. Roughly 
speaking, this is because, for every member of the statistical 
ensemble, there is another member, spatially displaced by an 

arbitrary amount. 

It is useful to compare the result ( 3.20| ) with the result nor- 
mally quoted for the spectrum rig„(/) produced in inflation 
(which is sometimes derived by assuming stationarity) [|[]: 
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(3.21) 



It is easy to see that in any practical experiment, o ne can not 
discriminate between these two possibilities: Eqs. (3.20) and 
(3.21 ). The reason is simple. In a practical experiment lasting 
(say) one year, the smallest frequency resolution df that can be 
attained is of order df w 1/year w 10^* Hz. This means that 
the energy density ( 3.20[ ) is resolvable in frequency bins not 
smaller than df. Hence the outcome of a (noise free) experi- 
ment would be a measure of pg^u averaged over a range of k 
for which dk — Andf /r]QHo. Even over this tiny range of fre- 
quency df, the argument of the squared cosine passes through 
a range d{kri) = Andf/HQ ^ 10^". The consequence is that 
in any practical experiment, the cos^{k'q) is averaged over ap- 
proximately 10^" cycles, after which it is indistinguishable 
from 1/2. Thus the stationary and non-stationary cases can 
not in practice be distinguished: the inflationary prediction of 
non-stationarity can not be falsified. 

Let us now return to the main line of reasoning and continue 
the calculation of the two-point correlation function and see if 
a two-point correlation experiment may be able to distinguish 



between the stationary and non-stationary backgrounds. For 
this purpose it is helpful to compare the inflationary model to 
a fictitious model Universe in which the particles were not cre- 
ated in perfectly correlated pairs but were instead formed by 
a stationary random process which did not correlate particles 
moving in opposite directions. In this case the wavefunctions 
have a time-dependent phase and are given by|TT| 



mv,k)] 



1 



Stationary 



V2 V Stt pp 



2r 11 (3 20) The se "m odes" should be compared with the ones given in 



Eq. (3.17). They lead to same average energy density in grav- 
itational waves as in the squeezed-state case. 

The two-detector correlation function C{t, t') for the infla- 
tionary model can now be derived simpl y by substituting the 
two-point function of the field operator (3.19) into the corre- 
lation function ( |3.4| ): 



5 
xrA r df -f{f)f^c^[r^{t)Mf)W[i(t')Mf)] 



(3.23) 



where the overlap reduction function 7(/) is a real function 
determined entirely by the relative separation and orientation 
of the two detector sites and is defined by Eq. (3.30) of refer- 
ence [0| (this function was originally defined and computed in 
reference p9[|). The overlap function 7(/) for the two LIGO 
sites is shown in Fig. Q The function 7(/) is unity for coinci- 
dent and co-aligned detectors. 
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** Many of the "standard" treatments of gravitational wave pro- 
duction by slow-roll inflation would imply that the energy-density in 
gravitational waves is zero. This is a misapplication of the standard 
consistency relation between the scalar and tensor amplitudes, be- 
cause in this non-tilted model the standard treatments imply that the 
energy-density in scalar perturbations is infinite. 



^^These are not normalized modes, since they do not correspond 
to vacuum fluctuations. Instead, they are modes containing the same 
energy density as in the inflationary case, but without the correlations 
between the oppositely-moving quanta. 



11 



FIG. 4. The overlap reduction function 7(/) for 
the LIGO-Hanford and LIGO-Livingston sites. The left graph has 
a linear frequency scale, the right graph has a logarithmic scale. 



Inserting the mode functions ( 3.17 ) yields 



iPds 



C{t.t')^T^i^ + Ze,r'!^Hi 






C. Are the non-stationary terms observable? 

It is now easy to answer the original question: could a coiTe- 
lation experiment carried out with two interferometric gravita- 
tional wave detectors distinguish between the non-stationary 
squeezed-state stochastic background produced by inflation, 
and a stationary background with the same average energy 
density? The answer is "no". The reason is the same as that 



df l{f)f ^ [ cos 27r/(t - t') + cos 27r/(t + t' + 4To)] , (3.24^iven following Eq. ^2^. Suppose that there were no signifi- 
cant detector noise to contend with, and that we were only try- 
ing to distin guish betw een the two possible correlation func- 
tions (3.24) and (3.25). In an experiment of realistic length 



where we have assumed that t and t' are in the present 
epoch, replacing k{rj — rf) by 27r/(i — t') and k{ri + -q') 
by 27r/(t + t' + 4To). As mentioned earlier, the number of 
cycles of a wave during the infinitesimal time interval dt at 
time t is dN — f{t)dt — kdri/2TT. Hence, provided that t] 
and T]' are not too far apart (in the cosmological sense) then 
one has k{rj — rj') « 27r/(i — t'). This holds provided that 
|i — t'l << Hq^. Note that, while it is tempting to replace k-q 
with 2-Kft, it is incorrect. In fact, provided that the inflation- 
ary phase is not too long, one has (today) k-q « 2TTf{t + 2Tq), 
where Tq = 2/3Hq is the present cosmological age. The "ad- 
ditional cycles" arise because as one goes towards the past, 
the frequency of the wave increases due to blueshifting, so, 
for example k{-q + rj') k 211 fit + t' + 4T^). 

In the stationary case, which lacks correlation between the 
amplitudes of the opposite-momentum modes, one would find 
that the correlation function is given by 



C{t,t') 



r(l + ^e,) 



~lPds_TT2 



Ihi: ' pp 

rf/7(/)r'cos27r/(i-i'). 



(3.25) 



This is identical to the inflationary case, except that the non- 
stationary term depending upon t + t' is absent. These expres- 
sions are valid provided that \t — TqI and \t' — To| are both 
very small compared to Tq. 

In a more complicated model of the early Universe, where 
the energy-density during the inflationary epoch was not ex- 
actly constant as here, but was instead a slowly varying func- 
tion of time, one would obtain an almost-identical result. The 
only difference is that one would find an extra, slowly-varying 
power law factor of /" in r2gw(/), and in the integrand of 
C{t,t'), where a is the so-called "tilt" of the spectrum. This 
slowly-varying factor would have no effect on our arguments 
or conclusions. 



**The apparent divergence of this integral as / — * is due to the 
approximation made in its derivation that the frequency / is between 
~ 1 mHz and ~ 1 kHz. The exact expression is free of infra-red 
divergences, but would give the same function of t and t' for any 
practical experiment, since the measured correlation function is given 
by the integral restricted to the bandpass of the detector. 



(say, one year) the smallest range of frequencies that would 
be observable is a bandwidth df — 10^^ Hz. Consider now 
the extra integral term / df ^{f)f~^ cos2t: f{t + t' + ATq) 
which distinguishes the two cases. The time t + t' + ATq that 
appears in this integral is approximately six times the total age 
of the Universe, in other words, about 10^'* seconds. Thus, 
even if the range of integration is restricted to 10^^ Hz, the 
cosine factor appearing in this integral undergoes more than 
10^ cycles. Since all of the other factors in the integral are 
smoothly varying over the range of integration, the resulting 
integral vanishes, in comparison with the integral containing 
the stationary contribution cos 2t: f{t — i') lB 

The same conclusion is reached when one uses the results 
and notations of Ref. [pi]]. The strong dependence on fre- 
quency of the integrand appearing in Eq. (14) of Ref. | pl] ] is 
due to the cos^[27ri/(f — f^)] factor, and not on the slowly de- 
caying t/2'^+i factor [The 2'Ki't of Ref. [ pi] ] corresponds to 
h] in our notation and 2TTvti, to fcryi]. We have already shown 
that 2TTvti, = krji << 1. However in the other term t is the 
cosmological age of the Universe, and 27ri/i — krjo > 10^^ 
is (27r times) the number of cycles that the wave has under- 
gone since the beginning of the Universe. Thus one cannot 
approximate the integral by taking its value at the lower limit 
of frequency, which is precisely the approximation that ap- 
pears to make the non-stationary nature of the process visible 
even in year-long experiments in Ref. [111]]. 

Finally we note that even if it were magically possible to 
observe the gravitational wave coiTelation over a sufficiently 
naiTow bandwidth for the second integral to contribute sig- 
nificantly, the best that one could possibly do with optimal 
filtering is to add a factor of 2 to the total correlation. This 
means that, even if the optimal filtering strategy for this non- 
stationary signal could be implemented, the most one would 
add to the total signal would be a factor of two. The total noise 
would remain the same. Hence, even if there were no prob- 
lems related to the short observation time (compared to the 
age of the Universe) it would not be possible to claim that the 
signal-to-noise ratio grows faster as a function of integration 
time for the squeezed background than for the standard one. 



^^ Note that in optimal filtering schemes to search for a stochastic 
background, one is mostly concerned with times t and t' for which 

\t-t'\ < 50 ms. 
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APPENDIX A: GRAVITATIONAL WAVE PREDICTIONS OF 
INFLATIONARY MODELS 



In this Appendix we review the predictions [ |12| , |14{ | of in- 
flationary models for the statistical properties of relic gravita- 
tional waves. We write the spacetime metric as 

ds^ = a{rif {^dff + [6ab + hab{jh x)] dx'^dx''] , (Al) 

where a, b run over spatial indices, 77 is conformal time, and 
0(77) is the scale factor For the gravitational wave modes rel- 
evant to ground and space based detectors, it is a good ap- 
proximation to take a{r]) to be constant at sufficiently late and 
sufficiently early times: 



0(77) 






(A2) 



for some rji and rjf. For simplicity and without loss of gener- 
ality we will take aj = 1. We define «&, (3k to be the Bogoli- 
ubov coefficients for the differential equation 



^^"{v) 



P 



a 
a 



Kv) = 0. 



(A3) 



In other words, ak and (3k are such that if ii{ri) — exp(ikri) 
for 7] < rji, then ii{rj) — ak e'xp(ikr]) + (3k exp(— i/cry) for 

77 > 77/. 

The metric perturbation hab for rj > 'qf can be expanded as 



habiV^^) 






{sA,M)-^B.m{fr) - -1^1 \Pkf Sab 5{f - /') 5\n, m) 
+ -\k\ alPk Sab 5{I + I') ^'(n, -m), (A5) 

TT 

where k = 27r/ and (5^(n, m) is the delta function on the unit 
sphere. If one drops the second term in Eq. ( k5 ), one obtains 



a stationary, Gaussian stochastic background with 



' ^gwiJ j 



fc4 



TT^Pc 



l/3fe 



(A6) 



Here rigw(/) is the usual energy density per logarithmic fre- 
quency in units of the closure energy density pc- The second 
term in Eq. (A5) encapsulates the non-stationarity. 

The stochastic background will be dominated by modes for 
which the number of quanta created per mode |/3fcp is large 
compared to unity, which for typical inflation models means 
all modes with frequencies / in the range 10^^® Hz ^ f (S 
MHz. For such modes the ratio of the coefficients of the first 



\ak\ 



l/5fc 



and second terms in Eq. ( A5 ) is very nearly a pure phase, since 



1. So we can write 



/ t -e\ t f'\*\ Pc"gw(,/ Jr 

(sA,n(/)sB,m(/ ) > - 4^2|J|3 ^^^ 



<5(/-/V^(n,m) 



(A7) 



27ri/nx A.n 
^ab 



where as before k = 27r/ and the phase x(fc) is given by 
exp[ix(fc)] — Cikl Pk- Hence, each inflationary model is char- 
acterized, in the large squeezing limit, by two functions of 
frequency: the spectrum 51gw(/) and the phase x(fc). Note 
that under changes 7/ — > 77 — At; in the origin of conformal 
time, x{^) transforms as x{^) ~* xi^) + "^kArj. Also the 
phase x{k) satisfies x(^^) ^ ^x{^^)■ 

In Sec. pu above we show that in a specific inflationary 
model, x(k) ~ 2krii, where 771 is a specific value of con- 
formal time around the inflationary epoch, and that fc77i ^ 1 
for modes that are relevant for ground and space based detec- 
tors. Thus we may take x{k) ~ 0. This conclusion is valid for 
all inflationary models. If one uses the method of calculation 
of Ref. [ P2| ] to approximately evaluate ak and Pk, one finds 
that, up to additive corrections of order unity, 

x{k)^2k7jk,E (AS) 

where rjk^E is the conformal time at which modes with 
wavenumber k re-enter the horizon, at which point 



(A4) 



k — a! I a. 



(A9) 



Here J d^iln denotes the integral over solid angles 
parametrized by the unit vector n, A runs over the two polar- 
ization components, and the tensors e"^" are the usual trans- 
verse traceless polarization tensors, normalized according to 
gacgbd^^^n^^^n _ 25j^g, We Specialize to a circular polar- 
ization basis for which e'^'^" = (e'^") . The quantities 
SA,n{f) are Gaussian random processes with SA.n(—f) = 
SA,n(/)* and whose statistical properties are given by [ po| , pi| ] 



For the relevant modes which re-enter during the radiation 
dominated era, 77^^; is independent of the details of the in- 
flationary dynamics, and from Eq. (A9) is given by rjk^E = 
?7o + 1/fc, where 770 is the conformal time such that a{rf) ex 
rj — rjQ during radiation domination (i.e. the extrapolated 
zero crossing of a{rj)). Thus from Eq. (A8) we find that 
x(fc) « 2 + 2^770. We can neglect the constant first term, and 
we can choose the origin of conformal time so that 770 ~ 0, 
thus giving x(^) = 0. This conclusion is used in Sec. || above, 
where we use a discretized version of Eq. ( A7) with x(fc) set 
to zero. 
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APPENDIX B: DETECTABILITY OF SQUEEZING IN A 
SIMPLE MODEL 

In this Appendix, we analyze the detectabihty of the non- 
stationarity in the context of the simple model of Sec. ||. Our 
starting points are Eq. (2.19), which descri bes th e statistical 
properties of the Universe modes, and Eqs. (2.22) and ( [2.23 ), 
which describe the relationship between the Universe modes 
and the measured modes. 



We can summarize the information in Eq. ( |2.19D in terms of 
a characteristic function; 



exp 



where 



N-l 



Y,{vfhj+v*h*) 



j=o 



exp 



-M^j) 



(Bl) 



N-l 



3=0 



(B2) 



and vi,. . . ,vn are arbitrary complex numbers. As before, 
e = is the stationary case and e = 1 is the squeezed case. 
We can derive the corresponding characteristic function for 
the measured modes Ha by using Eqs. (2.22) and ( |Bl[ ), which 
yields 



exp 



Af-l 



J2 {^aHa + s\Hl) 



A=0 



exp 



1 



Q{v, 



(B3) 



where now si, . . . , sm are arbitrary complex numbers and 
where on the RHS Vj is given by 



WatSa- 



(B4) 



Here and below it is assumed that repeated lower case indices 
j, fc, ... are summed over 0, 1, . . . , A^ — 1, and upper case in- 
dices A, B, ... are summed over 0, 1, ... M — 1 (see Sec. II C 



above). Now combining Eqs. (B2) and (84) yields 



6 = 'iTABSAS*B + 2e [xabsasb + XAss^Sg] 
where 



Bk 



is a Hermitian matrix, and 

XAB = E ^fe WAkWBk 



(B5) 



(B6) 



(B7) 



is a symmetric matrix. 

The difference between the stationary and squeezed cases is 
due to the matrix xab- It follows from the formulas (B3) and 
( |B5| ) that when xab = 0, the joint probability distribution for 
the variables Hi, ... , Hm is exactly the same in the squeezed 



and stationary cases; in particular the phase of each Ha is 
uniformly distributed over the circle. 

We now show that for a purely white stochastic background 
with (jj = constant, the matrix xab = and thus the squeez- 
ing has no effect on the statistical properties of the measured 
modes Ha, irrespective of the values of Ts and Tobs- First, 
using the formula (2.23) for the weights Waj we can obtain 
the formula for Tab '■ 



^AB — 



-y 



N-1 
^2m(AC-BD)/M S~^ J2 



X e 



-2tt%{C-D)j/N 



3=0 



which can be approximately evaluated to yield 

M 



Tai 






(B8) 



(B9) 



For the matrix xab we find from Eqs. (2.23) and (B7) the 
formula 



w-i 



XAB 



1 ^,2..(AC+..)/M^^. 



iV2 



2 -i-mjs/N 

3 ^ 



C,D 



3=0 



(BIO) 



Now in the case c r., = const of a white spectrum, the sum over 
j in the formula ( ^ICJ ) can be written as 



N-l 

E' 

3=0 



N 



1-a 
1-a ' 



where 



a = exp [-27^^(2js + C + D)/N] 



(BID 



(B12) 



We can assume that s > 0, which corresponds to T, > 0, 
since measurements today must have Ts of order the age of 
the Universe [P3|]. Also we can assume that s + 2M < N, 
which corresponds to Tg + 2robs < L, since we should take 
L ^ oo at the e nd of our calculation anyway. It then follows 
from Eqs. ^U\} and (|bT^) that xab = 0. 



Thus, the effect of the squeezing on the statistical proper- 
ties of the measured modes Ha depends on the spectrum aj 
of the stochastic background. Let us now turn to the case 
when the spectrum is colored. In this case it turns out that 
the matrix xab is small, and consequently the effects of the 
squeezing are small, whenever the observation starting time 
Ts is large compared to a correlation time that characterizes 
the spectrum of the stochastic background. In other words, 
the effect of the squeezing on the statistical properties of the 
stochastic background is time dependent; the effect is strong 
near t = 0, where the modes are all synchronized [|33|], but 
becomes weaker and weaker at later times. 

To see this, it is convenient to transform back to a contin- 
uum representation of the stochastic background. The spec- 
trum Sfii.f) of the background h{t) is related to the quantities 



aj by [M] 
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aj = ^N^AfShUAf), 
where A/ = l/(27rL). We define the quantity 



Chir) 



(B13) 



(B14) 



Then the real part of Ch [t) is the usual correlation function 
(h(t )h(t + r)) (in the stationary case). The discrete form of 
Eq. ( |b141 ) is 



ChijAr) 



iV2 



N-1 
fe=0 



'27Tl^ 



ShikAf), 



(B15) 



where A/ = l/(27rL) and Ar = 2ttL/N = LAt. Now 
by combining Eqs. ( ^13| ) and ( lB15|) we can evaluate the sum 
over j that appears in the formula([B10|) for xab'- 



E-l 



2^-^2Trt{2s+C+D)j/N 



-N^Ch[2T, + iC + D)AT]. 



(B16) 



The quantity (plq) will be small once Tg is much larger 
than the correlation time r» of Ch, that is, the time over which 
Chir) ~ Ch{Q)- This correlation time t* will be roughly the 
reciprocal of the shortest frequency scale over which Sh{f) 
has appreciable variation. The "high frequency" structure 
in the spectrum Sh{f) will presumably be dominated by the 
"breaks" in the spectrum 5/i(/) corresponding to transitions 
from one cosmological epoch to another (eg inflation to radia- 
tion domination). However, the physical time for such transi- 
tions to take place cannot be shorter than a local Hubble time. 
Hence, the corresponding correlation time t* will always be 
much shorter than the present day age of the Universe To, and 
consequently always much smaller than Ts for realistic mea- 
surements. Hence the matrix xab and consequently the effect 
of the squeezing on observational data will be very small. 

Turn, now, to the question of how accurately the ma- 
trix XAB can be measured, i.e., to the fundamental limita- 
tions on measurement accuracy imposed by "cosmic vari- 
ance". To address this question, we apply the analysis of Ap- 
pendix y to deduce the conditions under which the squeez- 
ing is detectable. Note that this analysis allows for the pos- 
sibility of combining the measurements of all the different 
mode amplitudes. Lets define the real random vector x — 
{KeHi , . . . , RcHm, Im_ffi, . . . , Im_ffj\/). Then, from Eqs. 
(B3) and (B5), the analysis of Appendix ^ appHes directly 
with 



Si 



Rer 

-imr 



imr 

Rer 



and 



Si 



Rex -Imx 
— Imx —Rex 



(B17) 



(B18) 



It then follows from Eq. (C9) that the difference between the 
squeezed and stationary cases is only detectable in the regime 
where the quantity 



A = tr 



(x-r-i) (x-r^ 



(B19) 



is large compared to unity. We can appro ximat ely evaluate 
the quan tity A by substituting Eqs. ( P9| ), ( |B10| ), and ( |B16| ) 
into Eq. (B19). This yields, in the limit L — > cx) and At -^ 0, 






(B20) 



A 



We now analyze the implicatio ns of the final result (B2C). 
First, it follo ws from the formula ( ^20| ) together with the defi- 
nition (B14) that in the regime Tobs ^ Tg, we always have 
< 1, irrespective of the nature of the spectrum Sh{f), 



which proves the claims made in the body of the paper. In 
the other regime where Tobs ^ Ts, we can see that A <C 1 
once Ts is much larger than the correlation time of t* of Ch- 
As argued above, t* will be much smaller than the present day 
age of the Universe Tq. Hence, experiments today would have 
difficulty detecting the non-stationarity even if they last a time 
-To. 



APPENDIX C: DISTINGUISHING BETWEEN TWO 

DIFFERENT DISTRIBUTIONS OF A GAUSSIAN RANDOM 

VECTOR 

In this Appendix we a ddress the following statistical issue, 
which arises in Sec. II C and in Appendix above. Suppose 
that X = {xi , . . . , xjv) is a zero-mean, Gaussian random vec- 
tor which satisfies either 



(XiXj) — 2^1 



(case 1) or 



iXiXj 



^2 ij 



(CI) 



(C2) 



(case 2). Thus, there are two possible variance-covariance ma- 
trices. Si and S2- In our application we will take case 2 to 
correspond to a squeezed stochastic background, and case 1 
to a stationary stochastic background. Suppose now that we 
have one measurement of x. How well can we distinguish 
between the two possibilities? 

We now show that, when Si — S2 <C Si, and when case 2 
actually applies, the two cases can be distinguished with high 
probability only in the regime where the quantity 



1 

-tr 

2 



(S2 • sr^ - 1) 



(C3) 



is large compared to unity, and not when this quantity is of 
order unity. 

It is easiest to address the question using the Bayesian ap- 
proach. Let the experimenter's prior probability for case 2 be 
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P2- Then after the measurement her probability for case 2 will 
be revised to p'2, where 



P2 



P2 



where 



^A _ P(x|2) 
p(x|l) 



(C4) 



(C5) 



and p(x| 1) is the probability that x is observed assuming case 
1, etc. It is clear that the difference between cases 1 and 2 is 
detectable in the regime A ^ 1, but not in the regime A ^^ 1. 
Inserting Gaussian probability distributions into Eq. (C5) 
we find that 

A = ix^- [Sj-i-Sa^i] •x+-lndet(Ei • ^2"^). (C6) 

Lets now assume that case 2 actually applies, so that the ex- 
pected value of XiXj is given by Eq. (C2). Then, typical values 
of A will be close to the expected value of the A which is 

(A) = -itr [1 - S2 • Sr'] - ^lndct(S2 • Sr')- (C7) 

Now suppose that the eigenvalues of S2 • Sj^^ are 1 + Xj for 
l<j<N. This yields that 



(A)-^E^-Mi + A,)] 



(C8) 



j=i 



If each I Aj I is small compared to unity, then to a good approx- 
imation we have 



= -tr I" 
2 L 



(C9) 



In our application we use the formula ( |C9| ) which should be a 
good approximation to the exact formula (Zl). 
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